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Abstract 

We apply a confining qq potential to charmonium and open charm states in order to model 
wave functions and to begin studying structure. Results (in momentum space) provide form- 
factor input to a four-flavor effective chiral Lagrangian which models dynamics of charmonium in 
hot hadronic matter. Estimates are made for J/tp dissociation cross sections and rates within a 
fireball. Our study attempts to improve on previous comover suppression models since it includes 
gauge-invariant form-factor formalism constrained by quark-model phenomenology. 
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I. INTRODUCTION 



Signals of momentary quark and gluon deconfinement in high-energy heavy-ion reac- 
tions are studied today more aggressively than ever before owing to current experimental 
activities at the Relativistic Heavy- Ion Collider (RHIC)Jl|. Hard probes represent a par- 
ticular piece to the overall puzzle whose goal upon assembly is to fully understand the 
strongly interacting many-body dynamics of ultra-relativistic heavy-ion collisions, complete 
with definite evidence for quark-gluon plasma (QGP) formation and propagation. The hard 
probes provide complementary information to such softer probes as photon and dilepton 
production. Very briefly for charm, the idea of Matsui and Satz0 proposes that in events 
where QGP is formed, screening effects tend to break apart the cc states leaving a "gap" 
between observed charmonium and expected. It is very convenient to look for evidence of 
charmonium breakup by studying the mass distribution of muon pairs and trimming away 
the background of non-J/tp contributions. There have already been suggestions that such a 
comparison might suggest glimpses of QGP||. 

Meanwhile, several authors have begun to systematically assess comover absorption using 
a variety of different approaches in order to improve understanding of the background due 
to such effects as light meson plus charmonium interactions leading to breakup [|, |[ 0, ||. |9|. 



10, O. Breakup of this type could possibly be misinterpreted as subhadronic effects when 



indeed, it is ordinary hadronic many-body physics. The aim of this study is to further refine 
estimates of hadronic cross sections for breakup of J/ip by constraining form factors with 
quark-model phenomenology and confining potentials, and then to use the form factors in a 
four-flavor chiral Lagrangian. 

Our article is organized in the following way. We discuss in Sect. II the confining potential 
and resulting meson wave functions. They provide information on the hadronic form factors 
to be later used in charmonium dynamics. Sect. [Hj includes a brief summary of the four- 
flavor chiral Lagrangian used to model the dynamics of the hadronic matter constituents. 
It also discusses gauge-invariant implementation of finite hadron size effects, namely form 
factors. Results are presented and discussed in Sect. |TV| and finally, Sect. |V| summarizes and 
concludes. 



II. CONFINING POTENTIAL 



Accomplishments of quark-model studies include rather detailed comparisons of calcu- 
lated hadron spectra versus observed for a long list of light mesons. The models are con- 
strained as firmly as possible by confronting such spectroscopic details as masses and decay 
widths. We take a recent result from the literature |3J and extend it beyond the light mesons 
to include D's r] c , J/ip and ip' (no isospin dependence is included, so we do not mention \c 
states explicitly). The potential is 

V K + Xr p + A ^ 7rK ' exp ( — r V r o) - - q\ 
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where the range r$ of the hyperfine term is taken to be mass dependent 
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A set of parameters is chosen to give the usual Coulomb plus linear form for the central part 

m u = md = 0.315 GeV; m s = 0.577 GeV; m c = 1.836 GeV 

k = 0.5069; k' = 1.8609; A = 0.1653 GeV 2 ; p=l 

A = -0.8321 GeV; B = 0.2204; A = 1.6553 GeV 5 " 1 . (3) 

By numerically solving the Schrodinger equation bound state wave functions are obtained, 
from which meson masses and rms radii are readily computed. Results for a selected set 
of mesons are listed in table |. Much of the motivation to do subnuclear structure in this 
way is to use the information as form-factor input to an effective Lagrangian description 
of the strongly-interacting hadrons. Typically, one views the spatial density as the Fourier 
transform of a (momentum space) distribution — the form factor. As is usually assumed in 
field theoretic modeling of this type, and is consistent with quark counting rules, a monopole 
structure is used. We extrapolate from massless to massive probes, and use a monopole- 
charge-form-factor-inspired expression (only meant as an indicator rather than a consistent 
model calculation) for the cutoff or off-shellness parameter 



The monopole form factor to be discussed later uses A in attempts to describe three-point 



vertices where a meson of mass m and rms radius y (r 2 ) is forced to go off shell. In rough 
terms, the size of the interaction vertex is inversely proportional to the cutoff parameter. 
And again, it is the off shell particle in a three point vertex which governs the physics. 



III. FOUR-FLAVOR CHIRAL LAGRANGIAN 



In the absence of firm experimental constraints on a full set of mesonic interactions involv- 
ing strangeness and charm, effective theories, with a specified chirally symmetric interaction, 
are quite useful. Indeed, there has been a renewed interest in these approaches since full 
understanding and control of the nonperturbative effects confinement necessitates is still be- 
yond grasp. Chiral symmetry and current conservation represent minimum requirements for 
any effective hadronic field theory. We take such an approach here by extending the usual 
two-flavor chiral Lagrangian to not three, but a four-flavor set of fields. The strange quark 
mass being greater than up and down quark masses probably already brings about some 
limitations for the Lagrangian's usefulness, and the charmed quark mass certainly breaks 
the symmetry to a deeper extent. And yet, it is not unreasonable to relegate this breaking 
to the mass terms, and insist that the interaction remain fully symmetric. 

The full details, starting from the nonlinear sigma model, introducing vector and axial 
vector fields into gauge covariant derivatives and then subsequently gauging away all of the 
axial degrees of freedom, have been published elsewhere 0. We therefore include here the 

TABLE I: Masses, rms radii and form-factor cutoff parameter for a select set of mesons. 



Meson Mass (MeV) [Obs.] 


Mass (MeV) [Calc] 


V(r 2 ) (fm) 


A (GeV) 


7T 


138 


138 


0.59 


0.80 


K 


496 


490 


0.59 


0.96 


P 


769 


770 


0.92 


0.93 


K* 


894 


904 


0.81 


1.07 




1019 


1021 


0.70 


1.23 


a\ 


1230 


1208 


1.24 


1.29 


D 


1867 


1862 


0.61 


2.03 


D* 


2008 


2016 


0.70 


2.12 




2980 


3005 


0.37 


3.25 


J/i> 


3097 


3101 


0.40 


3.32 




3686 


3641 


0.79 


3.73 



interaction terms alone. They are written very compactly as 

C[ nt = ^Tr(p,[9>,0])-^Tr([0,p^ 2 )+^Tr(^Kp1) + ^Tr(K,p1 2 ) (5) 

where and are 4x4 matrices with entries containing pseudoscalar and vector fields, 
respectively. 

Since a strict chiral symmetry is respected, there remains just one (chiral) coupling con- 
stant to fix. We do so by making certain the rho meson is correctly described. The choice 
g = 4.3 gives V p = 151 MeV at the pole mass of 770 MeV, and gives decay widths for 
strangeness and charm excitations listed in Table ||. 

TABLE II: Model prediction for the K* and D* decay widths as compared to experiment, 
particle Chiral Lagrangian Experiment 



K(892)° 44.5 MeV 50.5 ± 0.6 MeV 

K(892) ± 44.5 MeV 49.8 ± 0.8 MeV 

D(2007)° 10.1 keV < 2.1 MeV, 90% CL 

D(2010) ± 21.1 keV 96 ± 4 ± 22 keVp] 



The very recent charged D* decay measurement coming in at 96 keV allows for a D*Dw 
coupling constant to be fully twice as large as the chiral symmetry proposes. This would 
increase the dissociation cross section be precisely a factor of two. For now, however, we 
adhere to the chiral symmetry prediction. 

The interactions identified in Eq. @ do not include anomalous processes of type vector- 
vector-pseudoscalar. We therefore extend the set of interactions by introducing 

C£ t = g vv<t> Tr (e^W&'Vfy) , (6) 

with coupling constants constrained individually via vector meson dominance. One of the 
channels now open with £ n is J/ip + ir — ► r) c + p, an important contributor. 

A. Form Factors 



Effective theories attempt to model composite objects which necessarily have finite ex- 
tent, and are therefore responsible for finite-sized interaction vertices. Three-point functions, 



representing full details of the interactions including loop effects to arbitrary order, are no- 
toriously difficult to handle consistently. One-boson-exchange models, which involve three- 
and possibly four-point vertices, approximate the full calculation. Since in such models at 
most one particle per vertex goes off shell, a Lorentz invariant form factor accounts for dress- 
ing the vertex. So, rather than attempt to calculate higher-order contributions, we assume a 
specific monopole form for the momentum dependent vertex coupling "constants" , and use 
the off-shell or cutoff parameters from Table ID. All three-point vertices are therefore given 
the following monopole 

where t here is the squared four-momentum of the off-shell particle. Notice that when 
t — > m 2 , then h — *■ 1, which is indeed how the pole coupling constants are all defined. 

Four-point functions are modified from their typical g^ u form to the most general linear 
combination of all possible lowest-order Lorentz invariant structures constructible out of the 
external momenta. Specifically, in the reaction J /ip + tt — > D + D*, the four-point vertex 
becomes 

V i* = AgT + B (p>Lp»+p%p^)+C (^p'+ptp^ + D (p%p» + pLp^) 

+E(p» P : + p^p^), (8) 

and then the expansion coefficients A— E are chosen so that the overall scattering amplitude 
is fully gauge invariant, d^M.^ = 0. The choice is not unique; but on the other hand, gauge 
invariance alone is not enough to uniquely fix an amplitude. It does however represent an 
absolute minimum requirement of any reliable model. 



IV. RESULTS 



With all the interactions identified in the model, all the vertices constrained as much as 
possible, a list of reactions involving light meson plus J /if) can be enumerated and calculated. 
We begin looking at tt, K, r], p, uj, K*, and so on. 



A. Cross Sections 



The required dynamical quantity for estimates of J/ if) production and possible suppression 
are the breakup cross sections for the individual reactions. There are too many specific 
channels studied here to completely discuss them all. Instead, we show the two leading 
contributors in Fig. [I]. The pion-induced reactions involve D + D* , D + D*, and r/ c + p final 
states and the rho-induced reactions involve D + D, D* + D* and r\ c + n final states. In a 
rough summary, the cross sections range from a few tenths to a few millibarns. 



B. Dissociation Rates 



Consider the reaction J/if> + b — > 1 + 2. The rate within a fireball for this to occur is the 
following. 



dTj/^ = — !— d b f f f h | M\ 2 h h^Y5\ Vj/ ^ + Pb -p 1 - pa) (9) 
d 3 pi d 3 p2 

X ■ 



(27r) 3 2£! (27r) 3 2£ : 



2 



where db is the degeneracy of species b, fb is the Bose-Einstein distribution, and f = 1 + f 
to account for the medium. We show in Fig. |2] the total dissociation rates at two fixed 
temperatures as functions of the J / ip momentum. 



C. Flow 



It seems fairly clear by looking at the experimental results from RHIC that significant flow 



is present in the reaction zone[T3|]. Comover suppression is not expected to be significantly 
affected if the heavy charm is comoving. We suppose here, that it is not. We look at 
the possibility that J /if) breakup rates could depend on the radial flow velocity. There is 
no reason to expect the J/ips to have thermal momentum distributions since the elastic 
cross sections are too small to allow complete thermalization0]. We therefore estimate the 
new dissociation rates assuming rapid radial expansion, but with J /if) given a fixed three- 
momentum in the rest frame of the fireball. We then average over all solid angles with equal 
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FIG. 1: Cross sections for pion- and rho-induced dissociation. For complete details on all final 
states see Ref. M. 
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FIG. 2: Dissociation rates in a fireball at fixed temperatures, 150 and 200 MeV as a function of 
J/ijj (nonthermal) momentum. Rates included tt, K, p and iT*-induced breakup. 



weight. The expression we use "with flow" is therefore 



dT 



wf 



' pj/v- 
4ir 



dT 



(10) 



where dTji^ from Eq. (^) is used but the equilibrium distribution for particle h is now 
feqipb ■ U), where £7 = 7 (1, v) is the four- velocity. We use |v| = 0.6. Rates in the presence 
of flow are reported in Fig. |3]. 

V. CONCLUSIONS 

We have estimated charmonium structure using a confining quark potential calibrated 
with a long list of hadronic states. The momentum space wave functions provide form- 
factor input to an effective four-flavor chiral Lagrangian describing charmonium dynamics 
in a strongly-interacting many-particle system. Flow was introduced, albeit in a rather 
simplistic way, and was shown to affect the results. 

J/i/j breakup cross sections of several tenths to possibly a few millibarns were found. 
Kinetic theory was used to benchmark the dissociation rates in a fireball. At high J/ip 
momentum (5 GeV/c) and high system temperature (200 MeV), a dissociation rate of 10 
MeV was found. With flow present in the system, the dissociation rate increased by a factor 
of roughly 3. Future studies will include a folding of a more realistic pt distribution for <//?/>. 
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FIG. 3: Dissociation rates in the presence of rapid radial flow |v| = 0.6 assuming the J ftp has fixed 
three-momentum randomly distributed. 



